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ABSTRACT 

We propose that UV radiation can stimulate the formation of planetesimals in externally-illuminated 
protoplanetary disks. We present a numerical model of disk evolution including vertical sedimentation 
and photo-evaporation by an external O or B star. As solid material grows and settles toward the disk 
midplane, the outer layers of the disk become dust depleted. When such a disk is exposed to UV radi- 
ation, heating drives photo-evaporative mass-loss from its surface, generating a dust-depleted outflow. 
The dust:gas surface density ratio in the disk interior grows until dust in the disk midplane becomes 
gravitationally unstable. Thus, UV radiation fields may induce the rapid formation of planetesimals 
in disks where sedimentation has occurred. 

Subject headings: solar system: formation - (stars:) planetary systems: protoplanetary disks - (stars:) 
planetary systems: formation 
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1. INTRODUCTION 

The majority of stars are born from giant molecular 
clouds in dense but transient OB associations such as 
Orion. Nearby O and B stars can photo-evaporate the 
outer parts of circumstellar disks surrounding low-mass 
stars on 10 5 — 10 6 yr timescales. Several dozen disks 
with diameters from 100 AU to 1200 AU have been 
directly imaged around young stars embedded in the 
Orion nebula with the Hubble Sp ace T elescop e (HST) 
ijO'Dell WondlTTmri iBallv et a,l ill 99*3 Excess 
near-infra red emission i n the 2 to 10 /jm wavelength re- 
gion ijSmith et alJl2004|) provides indirect evidence that 
more than 80% of the several thousand low-mass young 
stars in the Orion nebula cluster are surrounded by disks 
of gas and dust. Disks which are not directly seen in im- 
ages are either too small, are lost in the glare of their cen- 
tral stars, or are hidden inside externally ionized cocoons 
('proplyds') produced by their own photo-evaporation. 
About half of these disks are being photo-evaporated 
by ultraviolet (UV) radiation from nearby O and B 
stars in the Trapezium cluster. Mass-loss rates have 
been measured to be in the range 10 -7 — 1O~ 6 M yr -1 
l|Hennev fc OTMll IT999T) . implying the loss of a disk 
within 10 4 to 10 6 years. It has been suggested that 
the formation of planets in regions such as Orion may 
be difficult - planets must form either rapidly, or rarely 
l|Throon et al.ll2001t iHester et alJl2004jl . 

In this letter, we use a numerical model to study the 
loss of material from a photo-evaporating disk and to ex- 
amine the effects of sedimentation towards the disk mid- 
plane on the loss of the gas and dust components. Pre- 
vious models of disk evaporation llJohnstone et aflll998t 
iMatsuvama et al.ll2003t lAdams et al.ll2004j) have exam- 
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ined only the gas component, assuming that gas and dust 
are well-mixed throughout and that the loss of gas im- 
plies the loss of dust. In our model, the evolution of the 
gas and dust components are treated individually. We 
explore the consequences of grain growth and sedimenta- 
tion on planetesimal formation. Our model assumes that 
planetesimal formatio n begins by non-gravi tational stick- 
ing in collisions {e.g. Cho kshi et al]ll993|) . after which 
material sediments to the midplane where it can un- 
dergo g ravitational instability af t er sufficient gas is re - 
moved HOoldreicn fc War HI 11975 fwiiTdin fc Shul 12002^. 
The model presented here builds on our earlier work that 
did not in clude vertical sedimentation or gra vitational in- 
stability ljThroodl2007i iThroon et al.ll2001ft . 

2. DISK MODEL 

We have developed a numerical model which tracks 
the state of the disk's gas and dust components as 
they undergo vertical sedimentation and external photo- 
evaporation. We consider the effects of UV-induced 
photo-evaporation due to nearby massive stars, and we 
ignore the effects of self-irradiation by UV light from 
the low-mass stars embedded within the disks. Thus, 
our model is most applicable to the outer portions of 
disks, or to disks which are effectively shielded from self- 
irradiation. 

Most disks around low-mass stars in OB associations 
formed well before photo-evaporatio n began. Low-mas s 
stars in Orion have ages 10 5 — 10 6 yr (Hillenbrand 1997), 
while most appear to have been pho to-evaporating; for 
only a few 10 4 yr ijBallv et alJ l2000h. Therefore, disks 
can evolve quiescently as grains grow and settle for up 
to about 10 6 yr, after which photo-evaporation begins. 

Initial Disk Str ucture. The init ial surface density fol- 
lows the model of iHavasliH l)1981|) . Surface densities for 
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the gas and dust as a function of distance r are set by: 
E d = 30.1 (r/AU)- fc [g cm- 2 ] (1) 



E 9 = 1700 (r/AU)~ fe [g cm" 2 ] 



(2) 



implying a dustrgas mass ratio of 1:60. We assume a 
mass distribution exponent k — 3/2 and a disk size r^ = 
100 AU, yielding a disk mass 0.02 Mq, with half of this 
between 1 AU and 30 AU . In the vertical direction z, 
we use the lHavashill(T98l|) pressure-balance scale height 
H = 7.1 x 10 n (r/AU)- 5/4 cm, and calculate the gas 
density 
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where the factor of two accounts for the disk's top and 
bottom sides. The central star has 1 Mq. 

Grain growth and sedimentation. Dust grains in young 
disks are transported by small-scale turbulent eddies. 
Within these edd ies, grains collide, stick, and grow 
(jCuzzi et al.lll996j) : as they grow, they sediment toward 
the disk midplane. Numerous studies have established 
that most of the disk's dust mass will have grown via 
collisions to sizes of cm or beyond an d sediment ed to 
the m idplane within a few 10 4 — 10 5 yr. Wcidcnschillins 
l|1997|) found that by 7 x 10 4 yr, typical bodies at 30 AU 
had grown to meter sizes and settled to the midplane. 
iMizuno et al.l l)1988|) found that typical g rains at 40 AU 
grew to nearly cm-sizes after 2400 yr, and lThrood {2000) 
found that grains grew to cm sizes at 100 AU within 

10 5 yr. HST observations directly support grain growth, 
even at the 600 AU outer edge of Orion's largest known 
disk l)Throop et al.1 12001^1 . Assuming a settling time 
t s ~ 10 6 (r/Vim)- 1 yr (e.g.. lYoudin fc Shdl200l . atyp- 
ical disk will easily be settled to beyond 100 AU in the 

10 6 yr for which it can evolve before the onset of photo- 
evaporation. 

Vertical structure. Local turbulence will prevent grains 
from settling to an infinitely thin plane. The equilibrium 
vertical distribu tion of dust in a proto-planetary disk has 
been studied bv iSekival l)1998j) . In their model, the den- 
sity is determined by an equilibrium between grain set- 
tling and grain lofting via the Kelvin-Helmholtz 'sand- 
storm' instability. Because the gas disk has only a finite 
mass-loading capacity, the dust density increases rapidly 
toward the midplane. In some c ases, the solution for m id- 
plane density becomes singular. lYoudin & Shul l)2002fi in- 
terpreted this singularity as being indicative of the dust 
disk's susceptibility to a gravitational instability (GI). 
They did not model planetesimal formation from this in- 
sta bility, on l y the onset of the GI. They found that in 
the Havashjl l)1981[) disks, an increase in the dust:gas sur- 
face density of ~ 2 — 10 (varying with r) was sufficient 
to reach this instability, and it could be caused either by 
increasing E<j or decreasing E g . 

The Sckiya vertical profiles assume that the disks have 
no global turbulence. This condition has not been estab- 
lished in the Orion disks; indeed, many models describ- 
ing turbulence from thermal convection and magneto- 
hydrodynamic instabilities have been pu blished (e.g., 
iLin fc PapaloizoulllQSOT: iBalbus et al]ll996fl . Because of 
their proximity to luminous sources of radiation, how- 
ever, Orion's disks are expected to be relatively warm 



with temperatures of around 50 K. Thus, turbulence pro- 
duced by thermal convection from the disk mid-plane is 
expected to be suppressed compared to that in similar 
disks in dark clouds such as Taurus where the disks may 
be as cold as 10 K. However, the precise nature of any 
turbulence remains an open issue. 

Observations pr ovide indirect evidence that the Orion 
disks have settled. iJohnstone et al.l l)1998f) estimate from 
model fits that the disks' outer shells are characterized 
by a UV penetration depth Nd — 3 x 10 21 cm -2 - that 
is, the outer shell is depleted in dust by at least a factor 
of 3 from the classical No- This would be expected in 
a disk in which solids had partially sedimented to the 
midplane. 

Photo- evaporation. UV radiation heats the disk sur- 
face to a depth Nd where the radiation is predomi- 
nantly absorbed by dust. Wherever the sound-speed 
in the heated layer exceeds the local gravitational es- 
cape speed, mas s-loss occurs. Mass-lo ss occurs in two 
distinct regimes ( John stone et aljl!998fl . FUV radiation 
(912A < A < 2,000A) produces a warm 10 3 K neutral 
hydrogen outflow for radii 77 > GM/2cj where cj ~ 3 
km s _1 is the sound-speed, and the factor of roughly two 
accounts for the effects of pressure gradients within the 
outflow. For a 1 Mq star, 77 ~ 50 AU. For most disks, 
the resulting mass loss prevents ionizing EUV radiation 
(A < 912 A) fr om reaching the disk surface. But, as a 
disk loses mass and shrinks to a size less than 77, FUV- 
induced mass-loss declines and stops entirely. The ioniz- 
ing EUV radiation can then reach the disk surface, ionize 
its skin, and raise its sound-speed to cji ~ 10 km s . 
Mass loss resumes in the form of a fully ionized wind 
until a radius of about rjj ~ GM/2c 2 II ~5 AU 



We incorporate the mass loss rates of IJohnstone et aT] 
(1998) from their eqs. (5) and (23). We assume an out- 
flow wind speed vq = 3 km s , and a UV penetration 
depth N D = 3 x 10 21 cm -2 (i.e., their e = 3). We as- 
sume that small, micron-sized dust grains produced as 
collisional byproducts will remain well-mixed with the 
gas and therefore continue to provide efficient heating 
of the gas as they absorb UV radiation, even as larger 
grains settle to the midplane. The external star is as- 
sumed to have an output of 10 49 photons s _1 and be at 
a distance d = 0.1 pc, roughly corresponding to Orion's 

182-413 (H ST10) disk. 

Following I Johnstone et ail l)1998l) . we assume that the 
mass-loss rate per area is constant across the disk. We 
assume that mass is lost from the top down (decreasing 
z), and that the mass loss occurs symmetrically on both 
sides of the disk. If the disk material re-settles vertically 
on a timescale faster than the photo-evaporation time, 
then the disk will maintain a general symmetry in the 
z direction, regardless of asymmetry in the illumination 
source. This condition is met (for dust and gas respec- 
tively), if 



and 
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Using typical grain sizes a(r) due to coagulative 
growth, and settling times t se t (e.g. IMizuno et al .1119811 
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IWeidenschiliiii3IT997t iThroonl 1200(1 . we find that these 
conditions are satisfied throughout our simulation. 

Dust entrainment. When gas is lost to the UV-induced 
outflow, small grains are entrained and dragged along 
with the flow. We assume, by analogy with the gas's 
finite mass- loading fYoudin fc Shi J 120021) . that gas can 
entrain its own mass density of dust grains, but no 
more. The entrained dust volume-density in the outflow 
is therefore pd e = Min(p,i(r, z), p g (r, z)). In most cases, 
these densities are equal near roughly iJ/100, causing a 
preferential removal of gas above this location. 

Model structure. Our numerical model uses a 100x100 
2D grid to track the densities p g (r,z) and pd(r,z). We 
use logarithmically-spaced bins in the vertical and ra- 
dial directions, and assume azimuthal symmetry. The 
model uses a self-adjusting variable timestep; typical 
steps are 10 3 yr at the outer edge and 10 yr at the in- 
ner. The timesteps chosen are more than sufficient to 
assure numerical convergence. At each timestep, photo- 
evaporation and grain entrainment act on the disk as 
described above. Our simulation begins when the disk 
first becomes exposed to UV radiation from nearby mas- 
sive stars; it ends when the disk has been eroded in to 
r<i < ru, stopping all photo-evaporation. 

3. RESULTS 

Figure n snows the output from our model run. The 
disk is gradually eroded from the outside edge inward, 
first by FUV radiation and then by EUV. After 1.25 x 
10 5 yr, the gas disk has been removed entirely outward 
of 5 AU, and continued UV flux will remove no addi- 
tional gas inward of this distance. Dust throughout the 
disk is partially but not fully removed: between 5 and 
50 AU, the disk is left with roughly 4 Earth masses of 
dust, or 3% of its original dust mass. This dust has 
condensed to the midplane and is not entrained in the 
photo-e vaporative outflow. I t meets the instability cri- 
teria of lYoudin fc Shul l)2002|K indicating that km-scale 
planetesimals could spontaneously form. Beyond 50 AU 
the dust density is not high enough to meet the instabil- 
ity criteria. Between 1 and 5 AU, an additional 20 Earth 
masses of solid material remains; this material is unaf- 
fected by external photo-evaporation and thus although 
it has settled, does not become unstable because gas mo- 
tions inhibit collapse. 

We have performed a second run on a disk with all of 
the same characteristics, except that the gas and dust 
densities £ s and E<z were increased to three times their 
nominal values, giving a disk mass of 0.06M©. In this 
run, the general behavior was similar to that shown in 
our nominal run. However, between 5 and 50 AU nearly 
140 Earth masses of material remained, or 40% of the 
original dust mass. The entire remaining disk met the 
instability criteria. The highly nonlinear dependence be- 
tween the disk mass and the fraction of dust retained is 
due to the fact that nearly all the additional mass goes 
directly to the midplane. Initial results indicate that this 
trend continues with even higher dust masses. 

4. DISCUSSION 

If a disk has more than a few 10 5 years to evolve be- 
fore it is exposed to UV radiation, grains throughout the 
disk can grow sufficiently large to settle toward the disk 
mid-plane. A gradient in the vertical distribution of the 



dust:gas ratio is established, and the disk surface layers 
become dust-depleted while the disk mid-plane becomes 
dust-enriched. Because photo-evaporation preferentially 
removes the gas component, the settled dust subdisk 
remains and the gas no longer inhibits a gravitational 
instability. The subsequent formation of planetesimals 
by gravitational instability addresses a size scale (cm to 
km) that has been difficult for purely accretional growth 
models, because bodies in this size range have neither 
the self-gravity nor the internal strength to maintain in- 
tegrity after collisions of more than several m s~ x (e.g., 
Weid enschilling fc Cuzzlll993th 

Planetesimal formation by means of gravitational in- 
stability may be enhanced by any process that sufficiently 
concentrates dust. For example, lYoudin fc Shul (2002) 
presented another method for enhancing the dust:gas ra- 
tio - the inward radial drift of cm-to-meter-scale bodies 
in the nebula. They demonstrated that local instabili- 
ties could be created on timescales of 10 5 -10 6 yr in the 
1-100 AU region. The methods presented in both our 
work and theirs operate on roughly comparable distances 
and timescales. 

Our model is simplistic and ig n ores some impor- 
tant processes. iHollenbach et all ()1994[) found that 
the central star's own UV radiation can drive photo- 
evaporation. Because we have ignored such self- 
irradiation in our models, the results presented here ap- 
ply to the outer regions of disks (beyond several AU), 
or to disks shielded from illumination from their central 
stars. Self-illumination acts in a similar way to external 
illumination and will also preferentially remove gas, so we 
expect that the effect in the inner region will be similar to 
that demonstrated here. Additionally, the inner disk will 
be affected by viscous transport of material toward the 
central star. Our model does not include these processes; 
iMatsuvama et all l)2003l) found that self-irradiation and 
viscous transport together formed gaps near the inner 
gravitational radius. Such gaps may halt the inward ra- 
dial drift of solid material, further spurring the onset of 
instability at that location. 

OB associations have been acc used of being haz- 
ardous to planetar y formation (e.g. iThroop et afll200H 
iHester et alJ l200"4Tl . On the contrary, our results indi- 
cate that UV radiation fields in an OB association may 
actually stimulate a critical step in planet formation - 
the growth of kilometer-size planetesimals. If our So- 
lar System formed via the method described here, the 
initial disk mass must have been greater than 0.02M© 
because that disk leaves insufficient material outward 
of 5 AU from which to form the giant planet cores, 
assuming the stand ard core-and-envelope scenario of 
iPollack et al.l l)1994[) . These planets have cores totaling 
roughly 40 Earth masses; our simulation leaves only a 
total of 4 Earth masses. If the disk were only several 
times more massive, however, the nonlinear dependence 
between initial and final dust masses would imply that 
sufficient mass for the cores would remain. The for- 
mation of planetesimals in our model depends on the 
loss of most of the disk's gas envelope, after which there 
would be insufficient material from which to form the gi- 
ant planets' atmospheres; therefore, if giant planets form 
near OB associatio ns, they ma y have formed by more 
rapid methods (e.g., Boss 2003). Partial retention of the 
gas envelope enabling giant planet formation could occur 
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if photo-evaporation ceases due to the short lifetimes of 
the illuminating stars, or if the gravitational radius r/j 
were moved outward from 5 AU. 

Our model may underestimate the amount of solid 
material left in the disk, for two reasons. First, as 
grains continue to settle out, the dust:gas ratio in the 
outer shells where photo-evaporation occurs will de- 
crease, causing a corresponding decrease in the amount 
of dust entrained. Because iVrj is considered a constant 
in our model, we do not account for this decrease in the 
fractional dust-loss rate. Second, our dust entrainment 
criterion has no grain size dependence. In reality, as 
grains grow they will become more difficult to entrain, 
and our method will overestimate their loss rate. In both 



cases, more detailed modeling of the gas-grain interac- 
tion is necessary; we estimate that after these factors are 
included, the disks may retain up to several times more 
material than indicated here. A future paper will present 
a more complete model of grain growth and entrainment. 
Additional work remains to be done to examine the sta- 
bility and evolution of volatiles in these systems. 
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Fig. 1. — Evolution of our model circumstellar disk of mass 0.02 Mq. Before the model begins, dust grains have grown and settled 
toward the midplane via Sekiya's vertical distributions. Photo-evaporation fully removes gas inward to 5 AU (left), while leaving substantial 
amounts of dust in this same region (right). The dust:gas surface density ratio through much of this region increases to a level where the 
dust is unstable against gravitational collapse (hatched region, right). 



